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Abstract 



Quark confinement and the genesis of tlie constituent quark model are exam- 
ined in nonperturbative QCD in Coulomb gauge. We employ a self-consistent 



> 

QQ I method to construct a quasiparticle basis and to determine the quasiparti- 

t^^ ■ cle interaction. The results agree remarkably well with lattice computations. 

1^^ ■ They also illustrate the mechanism by which confinement and constituent 

^^ , quarks emerge, provide support for the Gribov-Zwanziger confinement sce- 

-». \ nario, clarify several perplexing issues in the constituent quark model, and 

■^ I permit the construction of an improved model of low energy QCD. 
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I. INTRODUCTION 

Two of the key issues facing QCD at low energy are a quantitative description of confine- 
ment and an understanding of the origins of the constituent quark model. In this paper we 
demonstrate how both issues may be resolved through a nonperturbative analysis of QCD in 
Coulomb gauge. This demonstration makes the physical origin of both effects clear, resolves 
several longstanding inconsistencies in the constituent quark model, significantly extends 
the quark model, and establishes a - perhaps surprising - relationship between confinement 
and the constituent quark model. 

Although lattice gauge computations are capable of answering many questions in strong 
QCD, it is clear that the development of reliable analytical continuum tools are a necessity 
for advancing the field [^. Continuum methods allow one to understand how QCD works 
from first principles, permit the development of intuition for phenomenological model build- 
ing, and address computationally challenging phenomena such as QCD at finite density, 
extrapolation to low quark masses, or the treatment of large hadronic systems. A variety of 
such continuum tools exist: chiral perturbation theory, effective heavy quark and low energy 
hadronic field theories, 4-dimensional Dyson-Schwinger methods, fixed gauge Hamiltonian 
QCD approaches, and QCD sum rule methods. In this paper we focus on Hamiltonian QCD 
in Coulomb gauge. 

Much progress has been made in understanding Coulomb gauge QCD since the seminal 
work of Schwinger 0, Khriplovich 0, and Christ and Lee [Q. In particular, the problem 
of the Gribov ambiguity has been studied and a resolution has been suggested . The 
ambiguity arises because of residual gauge freedom after the canonical Coulomb gauge fixing 
condition, V • A = is imposed in a nonabelian theory. In Ref. [Q Gribov noted that 
the multiple-gauge copy ambiguity may be resolved by insisting that the Faddeev-Popov 
operator (to be defined later) is positive. With the aid of a simple model, he then showed 
that this constraint implies the existence of a novel form for the gluon propagator and an 
enhancement in the Faddeev-Popov propagator at low momenta. Furthermore, these imply 
that an enhancement exists in the instantaneous Coulomb potential, thereby providing a 
plausible mechanism for confinement. In a series of recent papers [§|,|^, Zwanziger has 
brought the Gribov Coulomb gauge confinement scenario onto firm theoretical ground and 
has demonstrated that a complete definition of the Coulomb gauge may be achieved by 
restricting the gauge fields to the 'fundamental modular region' - defined as the set of gauge 
fields which form the absolute minima of a suitable functional. Furthermore, the constraint 
to the fundamental modular region may be imposed by introducing a horizon term through 
a Lagrange multiplier in the Hamiltonian. 

A key feature of Coulomb gauge is that the elimination of nondynamical degrees of 
freedom creates an instantaneous interaction. The QED analogue of this is the Coulomb 
potential; however, the nonabelian nature of QCD causes this instantaneous interaction to 
depend on the gauge field, making it intrinsically nonperturbative for large fields. The re- 
striction of the transverse gluon field to the fundamental modular region formally makes 
the Coulomb potential well defined. It also implies that the Faddeev-Popov (FP) opera- 
tor which enters in the Coulomb potential is positive definite P,01- A consequence of this 
is that one may employ the variational principle to build nonperturbative models of the 
QCD ground state. This is a crucial step with many phenomenological repercussions in the 



methodology we will be advocating. As we shall demonstrate, the Fock space which is built 
on our variational vacuum consists of quasiparticles - constituent quarks and gluons. These 
degrees of freedom obey dispersion relations with infrared divergences due to the long-range 
instantaneous Coulomb interaction of the bare partons with the mean field vacuum. This 
interaction makes colored objects infinitely heavy thus effectively removing them from the 
physical spectrum. However, color neutral states remain physical because the infrared singu- 
larities responsible for the large self-energies are canceled by infrared divergences responsible 
for the long-range forces between the constituents. 

Constructing a quasiparticle basis is a nontrivial step which requires a nonperturbative 
treatment of QCD and, more directly, the QCD vacuum. We will show that it is possible to 
construct such a basis in a self-consistent manner by coupling a specific variational ansatz 
for the vacuum with the instantaneous interaction between color charges. The end results 
are explicit expressions for the Wilson confinement interaction, the spectrum of the quasi- 
particles, and the structure of the QCD vacuum. The resulting Fock space and effective 
Hamiltonian provide an ideal starting point for the examination of the bound state problem 
in QCD and provide a direct link between QCD and the phenomenological constituent quark 
model. 

A simplified version of this program has been investigated by the authors and others 
before P-[T^. In several of these studies the nonabelian Coulomb interaction was replaced 
by an effective potential between color charges, leading to a relatively simple many-body 
Hamiltonian with two-body interactions between constituents. The phenomenology of this 



approach has proven quite successful. In Ref. [|T^] we have extended this simple approxi- 
mation and treated the Coulomb kernel in a self-consistent way by considering the effect 
of resummation of a class of ladder diagrams. These diagrams originate from dressing the 
bare as/jx — y| Coulomb potential with transverse gluons. As one may expect from the 
discussion above, the effect of summing these diagrams is an enhancement of the Coulomb 
potential at large distances. Self-consistency appears in the problem because the strength 
of this enhancement is determined by the spectral properties the transverse gluons in the 
quasiparticle vacuum. 

In this paper we build on these findings by constructing a fully self-consistent set of equa- 
tions which describe the gluon dispersion relation, the effective instantaneous interaction, 
and the structure of the quasiparticle vacuum. A detailed derivation is given in Sec. II. This 
section also contains a brief review of the QCD Hamiltonian in Coulomb gauge and a discus- 
sion of the Gribov ambiguity. We discuss the renormalization procedure and show how the 
various counterterms in the regularized Coulomb gauge Hamiltonian may be constrained by 
physical observables. The last portion of Sec. II describes the variational vacuum employed 
in our method. Section HI presents the solution to the coupled equations. We first discuss 
the details of the renormalization procedure and present an approximate analytical solution 
which demonstrates many of the features which emerge. This is followed by a full numerical 
solution and a discussion of the effects of higher order terms. A comparison of these results 
to lattice data is presented in Sec. IV. Section V discusses the implications of our results for 
the constituent quark model and phenomenology in general. This includes clarifying several 
open issues in the CQM and extending the CQM. A comparison to similar approaches and 
our conclusions are presented in Sec. VI. 



II. QUASIPARTICLE FOCK SPACE FOR COULOMB GAUGE QCD AND 

CONFINEMENT 

One of the advantages of Coulomb gauge is that all degrees of freedom are physical. 
This makes the QCD Hamiltonian close in spirit to quantum mechanical models of QCD, 
for example the constituent quark model. The intuition gained from several decades of quark 
model calculations may then be applied to the analysis of a complex and nonlinear quantum 
field theory. Additional advantages of Coulomb gauge are that Gauss's law is built into the 
Hamiltonian, the norm is positive definite, and no additional constraints need be imposed on 
Fock space. Furthermore, retardation effects are minimized for heavy quarks; thus this is a 
natural framework for studying nonrelativistic bound states, and in particular for identifying 
the physical mechanisms which drive relativistic corrections, e.g. the spin splittings in heavy 
quarkonia. Since chiral symmetry is dynamically broken this framework is also of relevance 
for light flavors once the constituent quarks are identified with the quasiparticle excitations. 

The confinement phenomenon in QCD has two complementary aspects: (1) there is a long 
range attractive potential between colored sources; (2) the gluons which mediate this force 
are absent from the spectrum of physical states. Thus the mechanism for confinement is not 
particularly transparent in covariant gauges. In Coulomb gauge, in contrast, these two 
aspects can comfortably co-exist: the long range force is represented by the instantaneous 
Coulomb interaction and is enhanced as q^ — >■ 0, while the physical (transverse) gluon 
propagator is suppressed - reflecting the absence of colored states in the physical spectrum. 

A. Coulomb Gauge Hamiltonian 

Since the Hamiltonian in Coulomb gauge may look unfamiliar to many readers we briefly 
illustrate the derivation of the classical Hamiltonian here. 
The chromoelectric field is given by 

E'^ = -A" - VA°" + gT^^A^^A"", (1) 

and satisfies Gauss's law, 

V ■ E" + gr^^A^ ■ W = gp^g. (2) 

Here p^ = '?/'^(A"/2)^ is the quark color charge density. These equations are simplified by 
introducing the covariant derivative in the adjoint representation, 

j^ab ^ ^ab^ ^ Z^T^^feA^ (3) 

where T^ are the adjoint representation generators, T^^ = if^"-^. Thus Eq. (0) becomes 

D'^^ ■ E^ = gp^^. (4) 

If the electric field is split into transverse and longitudinal pieces, E = E^^ — V0 then Eq. (|[) 
yields 

- (D-^ ■ V)0 = <7p^ (5) 



where p"" = Pg + Pg is the full color charge density, with p" = /"^'^Ej,^ ■ A^ being the color 
charge density of transverse gluons. The equation of motion for the longitudinal component 
of the electric field, 

V ■ E" = -V ■ D"''A°^ = -V^0^ (6) 

leads to a constraint for the 0-th component of the vector potential which can be formally 
solved. This yields 



and 
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P"- (8) 



Finally the time evolution of the vector potential is determined by the transverse chromo- 
electric field through 

n-^ = -E^, = A'^ - g{l - V-2VV-)r*'A°''A^ (9) 

After canonical quantization, the transverse field 11" becomes the momentum conjugate to 
the transverse vector potential. A". 

Passing from the Lagrangian to the Hamiltonian yields terms proportional to (V0)^ from 
the longitudinal components of the chromoelectric field in E^, terms proportional to gpqA^ 
from the quark gluon vertex, gijj'y'^A^"'X"'/2ilj and terms proportional to gU'^ ■ A.^A^'^f"''"^ from 
the Ejr ■ Atr pieces of E^. Combining all these contributions and substituting the expression 
for A^ from Eq. (^ results in the instantaneous nonabelian Coulomb interaction. 

He = \j (fx(fyp\^)Kai>{^,y-A)p\y), (10) 



where 



ir,,(x,y; A) = (x,a|-/-(-V2)-/-|y,6), (11) 



V D' V D 

and p"" is the full color charge density as derived above, 

A_" 
~2 

The transverse conjugate gluon momenta 11" satisfy 



p"(x) = p^(x) + p^(x) = r''^A^(x) • n'=(x) + ^t(x)-V;(x). (12) 



A"'^(x), n^-^(y)] = iS'^' U^ - ^^] 5(x - y) = t6^%{V)i,6{^ - y). (13) 

Following Lee [|l^, we use the notation (x, a| ■ ■ ■ |y, 6) to denote kernels of integral oper- 
ators, 

(x, a|D|y, b) = [r^V, + ^r^''A^(x)l 5\x - y). (14) 



In the abelian limit D ^V, K ^ -^(^(x, a|l/V^|y, 6) = g'^6''^/A'n:\:s. - y| and the QED 
Coulomb interaction is recovered. 

A rigorous derivation of the nonabelian, quantum Coulomb gauge Hamiltonian was given 
by Schwinger and Christ and Lee ^. Zwanziger has shown how to derive the Coulomb 
gauge Hamiltonian with a lattice regularization ||^. The quantum Hamiltonian may be de- 
rived by transforming the canonical A^ = Hamiltonian to Coulomb gauge. The ^4° = 
Hamiltonian corresponds to 'Cartesian' coordinates in a flat gauge manifold, the subsequent 
restriction to Coulomb gauge induces curvature in the gauge manifold and therefore intro- 
duces a nontrivial metric. Christ and Lee have shown that the measure associated with this 
metric is proportional to the Faddeev-Popov determinant 

J = det(V • D). (15) 

Furthermore, the Hamiltonian contains factors of JT which are analogous to the Laplace- 
Beltrami operator induced when one first quantizes in curvilinear coordinates. The Faddeev- 
Popov determinant may be removed from the measure by working with the modified Hamil- 
tonian 

H -^ J^/^HJ-^/\ (16) 

which is hermitian with respect to ($|^E') = /X'A$*(A)^(A). Thus the final form for the 
QCD Hamiltonian in Coulomb gauge is 

H = H^ + H, + H^g + He, (17) 

where 

Hg= fd^ilj^-ia-V + f3m)ij, (18) 



H, = ^jd^ {j-'/'nj ■ nj-'/' + B • b) , (19) 



Hqg 



-g f dxtp^cx-Aip, (20) 



and 



Hc = lJ d^dyJ-'/'p^{x)j'/'Kat{^, y; A)j'/'p\y)J-'/\ (21) 

In order to compare with the covariant Feynman rules and the canonical path integral for- 
malism, it is convenient to Weyl order the operators (we note that Weyl ordering is the 
operator ordering which corresponds to path integral quantization with midpoint discretiza- 



tion). This leads to the Schwinger-Christ-Lee terms, Vi and V2 [17]. Here we will keep the 



original ordering of Eqs. (p!7|-21) so that no explicit Vi and V2 terms are present. 



B. The Gribov Ambiguity 



As detailed by Zwanziger |T8[, not only is the Haniiltonian renormalizable in Coulomb 
gauge but the Gribov problem can also be resolved ^ . The essence of the Gribov problem is 
that the condition V ■ A = does not uniquely fix the gauge in non-Abelian gauge theories; 
in general there are many copies of gauge field configurations, all with the same divergence, 
which are related by gauge transformations. Alternatively, the canonical transformation 
to Coulomb gauge is not singular so long as det(V -0)^0. But Gribov has shown that 
large gauge configurations exist such that this condition does not hold. As the true physical 
configuration space of a gauge theory is the set of gauge potentials modulo local gauge 
transformations, one must select a single representative from each set of gauge-equivalent 
configurations. The resulting subset of independent field configurations is known as the 
fundamental modular region (FMR). 

A convenient characterization of the FMR is given by the "minimal" Coulomb gauge, 
obtained by minimizing a suitably chosen functional over gauge orbits. This functional is 
defined as 

F^[g] = TTl£x{A3Y, (22) 

where g{x.) is a gauge transformation and A^ = gAg'^—g'Vg''. A simple calculation show that 
fields in the FMR are transverse. Alternatively, Zwanziger has demonstrated that Gribov 
copies may be removed by imposing the constraint {G)/V = (called the horizon condition) 
and argued that in the infinite volume limit imposing the horizon condition enables one to 
remove the direct restriction on the fields. Here G is the 'horizon term' given by 

G = y'dxdyD™(x) ■ (xa|;^^|y6) ■ B'^y) + (iVf - 1)V. (23) 

In this paper we follow a third approach. Because the Faddeev- Popov operator is positive 
semi-definite for fields in the FMR, we expand it in a power series over field variables and 
evaluate matrix elements by integrating over all fields. This is justified as long as the 
expectation value of the Faddeev-Popov operator does not change sign. We discuss under 
what conditions this procedure is consistent with the horizon condition in Sec. IV.B. 



C. Regularization and Renormalization 

To properly define the Hamiltonian a cutoff must be introduced to regularize ultraviolet 
divergences. This can be done, for example, by point splitting products of fields in the 
Hamiltonian. A simpler regularization procedure, adopted here, is to smear the fields. The 
induced nonlocalities are removed as the cutoff is taken to infinity. Since in the numerical 
studies to follow we will be working with renormalized quantities only (which are cutoff 
independent), we will explicitly remove the regulator making details of the regularization 
irrelevant. 

Counterterms need to be added to the canonical Hamiltonian to ensure that a cutoff 
independent spectrum is produced. 



H -^ H{A) -^ H{A) + SH{A). 



(24) 



In this paper we concentrate on the pure glue sector with at most static quarks, and therefore 
we will ignore the part of the Hamiltonian involving momentum or spin of the quarks. In 
the gluon sector, the presence of the cutoff leads to a single relevant operator (an operator 
whose canonical dimension is less then four). Thus SH{A) contains a term 



SH{A) = A 



2 ^m 



(A) 



(ix 



A"( 



+ 



(25) 



where Zm{A) is a dimensionless constant, and the notation [■ ■ -Ja represents the effect of 
regularization. For all marginal dimension four operators present in the canonical Hamilto- 
nian there will be corresponding operators in 6H{A) and the combination of the two leads 
to a Hamiltonian in which canonical terms are multiplied by A-dependent renormalization 
constants. For example. 



(ix 



n"fx)2 



+ 5 / rfx 



n° 



Zn{A) / dx 



n"fx)2 



(26) 



The full regularized Hamiltonian with counterterms is then given by 



H 



Zn{A) 



2 

^x(A) 



(ix 



n"(x)^ 



+ 



Za{A) 



(ix 



B" 



+ A^ 



,Zm[A) 



A^U? 



dxdy p'^(x)ir„,(x,y;A)p''(y) 



+ ... 



(27) 



The ellipsis stands for higher order terms induced by expanding the modified conjugate mo- 
menta i7~^/^njr^/^ in terms of gauge potentials. The effect of these terms will be discussed 
in Sec. III.F. 

At this stage we should in principle allow for every composite operator of dimension 
d appearing in the Hamiltonian to be multiplied by a renormalization factor Z(A)A^~'^ 
with Z being dimensionless and also allow for the coupling constant to be A dependent 
g -^ Zg{A)g = g{A). For example, as discussed earlier, if the fields are in the FMR the 
Coulomb kernel may be expanded in a power series in gA, and the order n contribution 
would be proportional to 



Zn{A) 



-A'T^ ■ V 



(28) 



Here Zn{A) is the n-th order triple gluon vertex (two Coulomb and one transverse) renor- 
malization constant and g{A) is the renormalized coupling. As we will show in Sec. III.F 
such vertices are UV finite which implies Zi{A) = 1. The contribution from the Coulomb 
kernel to the Hamiltonian can therefore be written in terms of only two renormalization 
constants Zk{A) and Zm{A) (and imphcitly g{A)) as in Eq. (|27|). 

As mentioned above, the A dependence of all renormalization constants has to be ad- 
justed in such a way that H leads to a A-independent spectrum. This implies that the 
renormalization group equations may be determined nonperturbatively from the spectrum 
of H. Furthermore in order for this Hamiltonian to be consistent with QCD (in the chiral 
limit) all renormalization constants Zi{A) cannot depend on A in an arbitrary way, but 
instead should depend on the scale through the coupling g{A). The renormalization group 
equations will be discussed in Sec. III.C. 
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D. Vacuum Structure 



The eigenstates of the Hamiltonian can, in principle, be expanded in an arbitrarily cho- 
sen complete basis which spans Fock space. One choice would be to use the perturbative 
basis which diagonalizes the free Hamiltonian H{g = 0). However, one expects the descrip- 
tion of any hadronic bound state would be very complicated in this basis. Alternatively, 
the phenomenologically successful constituent quark model indicates that hadronic wave- 
functions may saturate quickly with only a few Fock space states provided these states 
are constructed from constituent (quasiparticle) quarks. This strongly suggests that a ba- 
sis which incorporates the effects of spontaneous chiral symmetry breaking would be more 
efficient for describing hadrons and their interactions. 

We expect a similar scenario to apply to the gluon sector. In a given hadronic state there 
is a large probability of finding a component with a large number of bare, massless transverse 
gluons, but the expansion of a hadronic state may be significantly simplified in a transformed 
Fock space which is constructed from quasiparticle (massive constituent) gluons. We follow 
this intuition by constructing a vacuum upon which the quasiparticle basis is built with a 
functional Gaussian ansatz 



19 



^o[A] = {A\lj) = exp 



dk 

(2^ 



■A'^(k)cj(A;)A"(-k) 



(29) 



It may be shown pO[ that this ansatz sums all diagrams with nonoverlapping divergences. 
Note that the perturbative vacuum is obtained when u = |k|. The trial function is obtained 
by minimizing the vacuum energy density 

6 



6uj 



{to\H\u) = 0. 



(30) 



The vacuum state obtained from this procedure is denoted |ci;). We refer to u as the gap 
function since it is also responsible for lifting the single particle gluon energy beyond its 
perturbative value (see Fig. 5 below). 

This procedure is formally equivalent to the Hartree-Fock-Bogoliubov approximation, 
therefore one may also determine uj with a suitably chosen canonical transformation. Per- 
turbative gluon creation and annihilation operators are introduced in the standard way. 



A^(x) 



n-^fx) 




e(k, A)a(k, A, c) + e*(k, A)a'^(-k, A, c) 
— e(k, A)a(k, A, c) — e*(k, A)a^(— k. A, c 



„ikx 



„ik-x 



(31) 



with the perturbative vacuum satisfying, a(k. A, c)|a;(/c) = /c) = 0. The canonical transfor- 
mation is determined by requiring that the vacuum ansatz satisfies a(k. A, c)\uj) = 0, where 
the quasiparticle operators a, a^ are related to the fields by 

dk 1 r 

e(k, A)a(k, A, c) + e*(k, A)a'^(-k, A, c) 



A^(x) 



(2vr)3 ^2uik) 



,ik-x 



^(x) = -i 



dk uj{k) 



(2vr) 



e(k, A)a;(k, A, c) — e*(k, A)a^(— k. A, c) 



^jkx 



(32) 



The condition that emerges for uj{k) from Eq. ( PDD is identical to the condition that there 
are no a^a^ or aa operators in the full Hamiltonian. 

E. Self- Consistent Gap Equations 

The form of the QCD Hamiltonian in Coulomb gauge induces a crucial complication in 
the evaluation of the ground state energy density. This is because the interaction potential 
itself depends on the choice of the vacuum: the kernel K (Eq. 11) depends on the vector fields 
which depend on the gap function (Eq. 32). Thus the gap function is actually determined 
by a set of coupled equations which describe the vacuum energy density and the interactions 
which are used to obtain this energy density. This subsection describes how these equations 
are obtained; the solution is presented in the next section. 

The first step is the evaluation of the Coulomb kernel, Eq. ([TT|) . This is greatly simplified 



with the aid of the Swift equation |21 



Kab{:si, y; A) 



/(A) 



dgiA) 



(x,a|;^— ^|x,6). 



(33) 



The subscript A refers to the regularization of fields operators in the Coulomb kernel. Thus 
one need only evaluate the Faddeev-Popov operator (7/V ■ D to obtain the full instantaneous 
Coulomb kernel. This can be done by expanding the Faddeev-Popov operator in powers of 
gK and taking the appropriate contractions of the gluon field. The expansion is justified as 
long as the fields are restricted to the fundamental modular region. In the infinite volume 
limit, this restriction is not expected to affect field contractions as long as the expectation 
value of the horizon term vanishes. Thus the following expressions may be used. 



ah 



(a;| [A'^(x)A^(y)] Jo.) = -2 



dk 



^T(k) 
w(k;A)* 



^ik-(x-y) 



5afe (^ 



(a;| [n'^(x)n''(y)J^ k) = ^ / ^^'^T(k)a;(A;; A)e 



ik-(x-y) 



o; 



A'^(x)n''(y)J Jc.) = -(cl [n'^(x)A''(y)J^ |c.) = ^ ^ 



ah 



dk^T^^e. 



jk-(x-y) 



(34) 



We have temporarily allowed for A-dependence in the gap function. This is discussed in 
more detail in Sec. III. A. 

The expansion of the Faddeev-Popov operator is given by 

(xa||^|y,6)A = I)W(x,y; A)5., + ^ |rfziD«'=^(x, y, z^; A),, : A^^-^^-l,) :a 



Cl,«l 



Jn)ci---c„ . 



+ ... Y. E y cizi---rfz„A\"T""(x>y>Zi,---Zn;A)„,:A'=^'*Hzi)---^'"""K):A+... 

C.\"'C.n i\ ■■■in 

(35) 

where :: stands for normal ordering with respect to |a;) and c„ and in refer to color and 
spatial components of the gluon field respectively. Here D^^^ stands for the expectation 
value (VEV) of the Faddeev-Popov operator in the ansatz vacuum. 
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.(0) 



Di;:^(x,y;A) = (c.|(xa| 



^(A) 



\y,b)\uj)A 



(36) 



V D 

An operator expansion of the Coulomb kernel may be defined in a similar manner 

Kati^, y; A) = ^^^(x, y; A)^,^ + J2 J ^Zi<'^''(x, y, Zi; A)ai, : A^^'^Hzi) ^A 

+ ••• E E /rfzi---rfz„irf;:).7;""(x,y,Zi,---z„;A)„, :A^^'^nzi)---^'"""W 

(37) 

The equation for the VEV of the FP operator is most easily expressed in terms of its 
Fourier transform which we write as 



(27r) d{P) ^^ dab 



rfxrfyDi?(x,y;A)e^'^-(''-^)e^^-^. 



(38) 



The amplitudes Z^*^"-* which multiply a product of n gluon fields can be written in terms of 
D^^' and a set of vertex functions, F*^"'. To do this we first define the Fourier transform of 
the £>(") via 

(27r)3<5(k - E q. - l)dffX'"^{k, qi, ■ ■ ■ , q„, 1; A)= - J d^dydz, ■ ■ ■ dz^e-^'^-^+^'-^^'^.t, ^--^ 



,(n)ci---c„ 



xDi:^!X (^'Zi'---'Z"'y;A) 



(39) 



Next we define the full transverse gluon-Coulomb vertex as F^(k, q, p). The Dyson 
equation for the full vertex is illustrated in Fig. 1 and is given by 



r^^Hk,I,k-I;A) 



n(k,q,k-q;A)^Z,(A)T^.^ + EEf (^^ 

X r^(k - 1, q, k - I - q; A)^^^i^^r^;(k - 1 - q, 1, k - q; A) 



c/(k - I; A) 
(k - 1)^ 



(40) 
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FIGURES 



nnn + 




I riyir^ 



FIG. 1. Dyson equation of the Coulomb-transverse gluon vertex. The thick line represents the 
full FP function (i(k). The hatched circle represents the full vertex F?. The gluon line is the gluonic 
quasiparticle. All external propagators are truncated. 

In the planar approximation higher order vertex functions F"^") satisfy the following Dyson 
equation, 



T(n);ci---c„ 



(k, qi, ■ ■ ■ , q„, k - ^ Qi) = ^i^.Z"''"{K qi, ■ ■ ■ , Qn, k - ^ q, 



j=i 



i=l 



dqo 1 



+ 5^7 (27r)3 2c.(go) 



r^?(k,qo,k-qo 



d{k- qo) ~(„).e^...c„ 
(k-qo 



,2 n---«n 



(k - qo, qi • ■ ■ q„, k - qo - XI Qi 



j=i 



X 



d{k - qo - EILi (Ir. 



(\c c V" ^ N2r^?(k-qo-Eq«'qo,k-Xq„ 

where we have introduced the following quantity: 
r&-'^"(k,q„...,q„,k-X:q,)=EE ^ ^"^^ 



(41) 



«1---«T1 



j = l 



Co «() 



(27r)3 2a;(go) 



F?«(k,qo,k-qo: 



xF-Kk-qo,qr,k-qo-qO ^-^^"^°"^^^: x 



(k-qo - qi) 



n-l 



r-;(k - qo - E ^i' qn, k - qo - E qi 



j=l 



1=1 



^(k-qo-Er=iqn) 

(k-qo - ELi qn)^ 



X 



r^(k - qo - E q«' qo> k - E q^ 



i=l 



i=l 



c^(k - qo) 
(k - qo)2 



(42) 







+ 
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FIG. 2. Dyson equation for F'^^. Symbols are as in Fig. 1. 



The equation for F*^^) is shown in Fig. 2. 

Finally, we are able to write the coefficients of the operator product expansion of the 
Faddeev-Popov operator as 



v(i) 



C (k, q, k - q) = d{k)mk, q, k - q)rf(q) 



(43) 



and 

4ilc2;ii,i2(k, qi, q2, k - qi - qa) = c/(k)F"i(k, qi, k - qi)c?(k - q^F'^^j^ _ ^^ ^^ k - qi - qa) 

+ dikjI^^^-^'^^Ck, qi, qs, k - qi - q2)rf(k - qi - ga), 

(44) 

and similarly for higher orders. Before renormalization, these amplitudes are functions of 
the cutoff. In the planar approximation the VEV of the Faddeev-Popov operator, d{k; A) = 
rf(°)(k,k; A) defined in Eq. (|8D satisfies. 



d{k;A) 



^(A) 



where 



^[^-]^?]^^^5.^l^/ 



l-g{A)I[d,uj] 
1 [^ dqi dq„ 



(45) 



n 



(27r)3 (27r)3 2cj(gi;A) X ■■■ X 2cu(g„;A) 



i(n)ci---c„/i , \^ „ \ ^ Z.^i=iqi;A) („)ci---o 



:nrX'"(k;q.,---,qn;k-5:q, 



j=i 



(k - ELi q.) 



nrx'"(k-Eq^;qi'---'q-k). 



1=1 



(46) 



The trace is taken over the implicit color indices of the vertex functions, F*^") = T^ , which 
also absorb the renormalization constants Zi{A) of Eq. (PB]). This equation is shown in 
Fig. 3. 



g 1 - 8 




FIG. 3. Dyson equation for the VEV of the Faddeev-Popov operator. See Fig. 1 for an 
explanation of the symbols. 
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We proceed to the evaluation of the Coulomb kernel. Following Swift [B^ we define 
/(k;A) via 



,/(k,A) 



{2nf6{P)d\K,A)^^^6,, ^ Zk{A) J d^dy{u\K,,{^,y)\u)^e^^<''-y^e'''' 



(47) 



From Eqs. (|33D and (|45|) it follows that the vacuum expectation value of the Coulomb kernel 
satisfies 



f{k,A) = ZK{A) + ^I[ZKd,cu]. 
dg 



(48) 



This comprises a linear integral equation which must be solved for / after having obtained 
d. We are finally in a position to evaluate the expectation value of the energy density from 
the full Hamiltonian, 



c- On \ Orn ~r c-i 



1 



C 



ViN! 



-{uj\H\uj) 



(49) 



where the three terms represent the kinetic energy (including the nonabelian portion of the 
B^ term), the mass counterterm, and the Coulomb potential respectively. In particular, 



Sn 



and 



1 f^ dq 
2J (2^ 



Zn(A)^(g;A) + Z^(A) 



uj{q;A) 



2,,.N, /-A dq dk (3-(k-q)2) 



+ ^'(A) 
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(27r)3(27r)3cj(g;A)cj(fc;A) 

(50) 






\z A' f""^ ^— 

2 "^ J i2nyu{q;A)- 



(51) 



The contribution from the Coulomb potential may be evaluated with the aid of the 
operator expansion in Eq. (pT]). Recall that the products of gluon fields in the operator 
expansion of the kernel are normal ordered with respect to the variational vacuum. Thus 
the maximum number of terms which contribute to the vacuum energy energy density is 
determined by the number of external fields present in the charge densities multiplying the 
kernel {i.e., four). The Coulomb vacuum energy density may be thus be written as 



Sc = S. 



(0) 

c 



p(2) , p(4) 



(52) 



The terms f^'-"^ correspond to the vacuum expectation values of K^^'' contracted with the 
fields from the charge densities. For the first term one gets. 



c(0) ^ Nc 
16 



^c 



/■A dq, dk /(k + q;A)rf^(k + q;A) ^ 

J (27r)3(27r)3 (k + q)2 I + I " ^J 



' uj{k;A) ^ uj{q;A) _ ^ 
uj{q;A) uj{k]A) 



(53) 



The higher order terms 8^ are of order d^^^'^\k] A). Since d plays the role of the running 
coupling (see Eq. (|^)), we expect these higher order terms to give finite corrections to 
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uj{k; A) which should be small in particular for large momenta because of the suppression 
from the running coupling. The effect of these higher order terms as well as vertex corrections 
(c.f. Eqs. (^) and (^iD ) and the FP determinant will be discussed in detail in Sec. III.F. 

Minimizing £(j^ with respect to u leads to two contributions - one from the explicit to de- 
pendence (cf. Eq. ( |5B| ) for Sq ) and the other from the implicit uj dependence arising through 
the kernel f(P. We refer to these contributions to the gap equation as 8^ '^ and 8^ re- 
spectively. The first of these is of order 0{S^^'^\k] K)) and the second is 0{S'^^'^\k]K)). 
Thus, for example 8^ should be combined with other order 0{d^{k]K)) contributions 
from 8(j '^ . Subsequent expressions for 8q "^ contain a factor of —2ixP' with respect to 
the derivatives of 8. For the moment we retain only the leading (9((i^(A;; A)) contributions 
from 8q'^ in the gap equation. Minimizing 8 with respect to oj leads to the following gap 
equation 

Z^(A).^(. A) ^ Zim + ^™(A)A^ + /(A)^ f ^^^^i^ + 

AT, M rfk /(k+q^AK(k+q^ ^ oA^^A^^u^^A) 

(54) 

This completes the derivation of the leading order gap equations. To summarize, these 
comprise Eq. P5| ) for the VEV of the FP operator dik,] A), Eq. (Q) for the Coulomb kernel 



f{k\ A), and Eq. ( ^if ) for the gap function ijj{k\ A). 



III. SOLUTION OF THE SELF-CONSISTENT GAP EQUATIONS 

Before continuing we shall briefly summarize our philosophy. The goal is to construct 
a quasiparticle Fock space which will provide a useful starting point for the evaluation of 
hadronic observables. Quasiparticle states are built on a variational vacuum and reflect the 
propagation of these degrees of freedom through a nontrivial background. Of course the 
full Hamiltonian still contains many-body terms which mix the free quasiparticle states; 
nevertheless, the quasiparticle Fock space is complete and at least in principle one should 
be able to diagonalize the full Hamiltonian in this basis. 

When dynamical quarks and gluons are considered, one would need to diagonalize the 
Hamiltonian in the full Fock space. In practice, however, such diagonalization is always 
performed in an appropriately selected subspace e.g. including only \QQ) or \QQQ) quasi- 
particle states. Such a truncation is better justified when the quasiparticles behave as 
constituent particles with average kinetic energies of several hundred MeV. Furthermore, as 
discussed earlier, the quasiparticle basis diagonalizes the one-body part of the Hamiltonian, 
thus at least at the level of Tamm-Dancoff truncation, the quasiparticle vacuum decouples 
from the hadronic spectrum. 

The required cut-off independence of the eigenvalues can be used to determine the A 
dependence of the various counterterms and couplings. At this stage, this implies that 
the Fock space itself should be cutoff-independent because, for example, the ground state 
energy of two static color sources is directly related to the expectation value of H in the 
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variational vacuum. We note that expanding the Fock space in which the Hamiltonian is 
being nonperturbatively diagonahzed will add new counterterms to the Hamiltonian which 
will modify the renormalization group equations. 

A. Vertex Truncation 

We start by examining the renormalization group structure which follows from the re- 
quirement that the gluonic Fock space itself is A independent. This implies that the Coulomb 
kernel, and hence /(/c; A) and d{k]A), should be A-independent. These conditions may be 
imposed through an appropriate choice of the cutoff dependence of the counterterms and 
coupling. 

Consider first renormalizing the FP operator of Eq. (^31). In this equation d is expressed in 
terms of the vertex functions F*^"^ and the gap function u. Since these can be independently 
renormalized using other renormalization parameters which do not explicitly show up in 
Eq. (^) {i.e. Zi{A), Zm(A) in Eqs. (^) and (pH])) we can replace them by their renormalized, 



A-independent versions, F^"^(---;A) -^ F^"''(---) and uj{k]A) — >■ uj{k). Thus the only A- 
dependent parameter available to enforce the cutoff independence of the FP operator is the 
coupling, g{A). 

To determine the consequences of this observation we examine the behavior of the vertices 
which appear in the equation for d (|45|). Asymptotic freedom implies that for momenta 
near the UV cutoff, the gap function and the renormalized vertex functions approach their 
corresponding free-field values. 



lim uj{k) = k + 0{g\A)) (55) 



and 



r^(k, q, k - q) ^ k^r + 0{g\A)). (56) 

For n > 1 

n 

rW (k, qi, ■ ■ ■ , q„, k - ^ q,) ^ 0{g-+\A)). (57) 

Similarly one expects that in this limit d{k) — > 0{g{A)). Thus the integral in Eq. |^ is 
logarithmically divergent as A ^ oo. This divergence is absorbed by the coupling g{A). It 
follows from Eq. (^) that F^ is given by an expression which is finite as A approaches infinity, 
thus there is no need for vertex renormalization and one can set Zi{A) = 1. Furthermore, 
the correction to the bare vertex T'^k is expected to be of the order 0{{g^)) where {g) refers 
to an UV and IR finite integral over the running coupling. This is due to the two Faddeev- 
Popov operators d{k) in Eq. (^Ol). Since d{k) is proportional to g{A = k) for large k, the 
renormalized FP operator can be associated with the running coupling: 

lim d{k) -^ g{A = k). (58) 

k^oo 



From the resummation implicit in Eq. (^5[), and consistent with asymptotic freedom, the 
large momentum behavior of d{k) will be logarithmically suppressed with k. Furthermore, 
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if d{k) is less singular than 1/k in the infrared limit then the integral on the right hand side 
of Eq. pop represents a finite, higher order (in the running QCD coupling) correction to the 
bare vertex. This is also true for the higher order irreducible vertices, T^"\ From Eq. (^iD 
it follows that these are 0{{g'^~^'^)). This important observation will be used to truncate the 
gap equations in the next subsection. 

B. The Truncated and Renormalized Gap Equations 

The considerations of the previous subsection may be used to truncate the general gap 
equations derived in Sec. II. This is necessary to make the equations tractable. The effect 
of neglected terms will be discussed in Sec. III.D . 

We start by ignoring the finite higher order corrections to the vertices and thus take 

n(k,q,k-q) = T'=F, (59) 

and 

r(") -^ 0. (60) 

The equation for the unrenormalized FP operator, Eq. (|i5|), becomes 
1 1 .. /"^ dq l-(k-q)2 ,^, 



dik;A) giA) ^J (27r)3 2^(g)(k - q 

One sees from this equation that in order for d{k; A) to be A independent, g{A) must obey 
the following renormalization group equation 

+ ^c / 77^-T^ TT-T-VT- -^diq - Ai) . (62) 



giA) gifi) V (27r)3 2.;(g)(q-M)2 

Thus Eq. (^) becomes 

' - ' ^c/,^^^y^^(q-k)+ivJ^jy^^.(q-M). 



d{k) d{fi) "7 (27r)3 2Lu(g)(q-k)2 ^" ' '' J (27r)3 2cu(g)(q - /^)2 

(63) 

Here the renormalized FP operator is written as d{k] A) -^ d{k). Eq. (^) implies that d{k) 
is independent of A (and the scale /i), and represents the once-subtracted form of Eq. (|6lD. 
The presence of g{A) in Eq. (^5|) shows that d{k) can only be determined up to an overall 
constant. Thus the equation for d{k) contains a single unknown, d{fi). 

The vertex truncations and Eqs. (^Sf) and (|63|) imply that the expectation value of the 
unrenormalized Coulomb kernel is given by 

nk,A) = Z.(A) + N.lj^, ^-J^^X/ f(.l - k)/(q ^ k; A). (64) 

The UV divergence from the integral is absorbed by Zk{A). Subtracting once yields 
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/(*) = /(.) + K I ^^igi-^.= (q - k)/(q - k) 

Here /(/i) is another external renormalization parameter. The renormahzation constant is 
given in terms of it by 

Z.W = /(.) - N.J^^^ZM^,^^ _ ^,y,q _ ^,. (66) 



We finally discuss renormalization of the gap equation, Eq. ( p^) . In general this equation 
can depend on the three renormalization constants, Z^(A), Zn(A), Zm{A.) and the renor- 
malized coupling, g{A). The coupling g{A) is already determined by Eq. ([6^). In the UV 
limit the integral on the right hand side of Eq. (Q) has in principle quadratic and loga- 
rithmic divergences. The logarithmic divergence is present if the kernel /(k — q)(i^(k — q) 
approaches a constant in the UV limit. There are, however, logarithmic corrections to both 
/ and d which follow from Eqs. (|63|) and ( |65D which actually protect the integral from the 
logarithmic divergence. Thus one can immediately set Za = -Z'n = 1 and absorb all possible 
remaining divergences (as A -^ oo) into Zm{A.). This leaves the quadratic divergence which 
is eliminated by a single subtraction, 

2 2 , 2, , , iVe / dk /(k + qK(k + q),^ , ,,^ ^^,^uj\k)-uj\q) 



c.(g) = g -/. +^(/^) + xy (^ ^^qr^p (l + (k-q.. ^(^^ 

k-^) ) — . (67) 



4 J (27r)3 (k + //)2 ' ^ ^' ' iu{k) 

The mass counterterm is given in terms of ui{fj,) by 

Nc /-A dk (3-(k-q)2 



ZUA)A' = u:\f.) - f^' ~ 9\A)-^ J 



(27r)3 uj{k) 



AT, ^A dk /(k + ^M^(k + ^),^ ^ ^^^^ ^^,^ uj\k)-u;\f,) _ ^gg^ 



4 y (27r)3 (k + /x)2 ^ ' ^ ^^ ^ aj(fc) 

Equations (|63|), (|65|), and (^) form the renormalized coupled gap equations which rep- 
resent the leading order vacuum and quasiparticle structure of QCD in Coulomb gauge. We 
proceed by examining the perturbative limit of these equations before turning to analytical 
and numerical solutions. Subsection III.F examines corrections to the gap equations due to 
truncation to the leading terms. 



C. Asymptotic Renormalization Group Equations 

We establish the relationship of the renormalized gap equations to standard perturbative 
QCD in this section. The renormalization group equation for the renormalized coupling, 
Eq. ( |62D implies that for large cutoffs 



JgjA) _ 8Nc 9'{A)d{A) 

dA ~ 3 (47r)2 ' ^''''^ 

and from Eqs. (|62D and (|63|) it follows that in the limit A — >^ oo 

We call the first coefficient in the expansion of the jS function /3q. The last equation 
implies that 

A ^ 'f. (71) 

Although it is tempting to compare this to the canonical perturbative expression of /5o = 
llA^c/3, this is misleading for two reasons. First the coupling defined here corresponds to 
the product of the VEV of a composite operator {i.e. the Faddeev- Popov operator) and 
the QCD coupling. Thus (3 will also refiect renormalization of the FP operator. We note 
that this is nevertheless a sensible definition for the coupling since it is this product which 
determines the strength of the various interactions involving Coulomb gluons. The second 
reason is that we sum loops which arise from the expectation value of the Hamiltonian 
and do not include those from iterating the Hamiltonian. Iteration of the Hamiltonian 
involves summing over intermediate states. This is fine in perturbation theory, but because 
of confinement can only be justified for color singlets so that summation should be restricted 
to hadronic intermediate states only. As discussed in Sec. H.D this may be achieved in bound 
state perturbation theory once the quasiparticle Fock space is specified. We will discuss the 
running coupling in more detail below. 

The expression for Zk given in Eq. ( |68D implies that the renormalization group equation 
for Zk{A) is given by 

Z^(A)_ rf^(A)/(A) 
^^A~ - ~^° (47r)2 ' ^^2) 

which in the limit A — > oo leads to 

. dZKJA) 8Nc9\A) 



Finally, Eq. (|68D yields 



^dZUAl_ 2^ .^. ^^(A) 



'^+ii!|^^°)+^-(^) 



(74) 



The first term is universal and refiects the quadratic divergence. The remainder relates 
to the UV behavior of the Coulomb kernel and the quartic-gluon vertex which are both 
determined by the running coupling g{A). 

As expected, all counterterms run as a function of a single renormalized parameter g{A), 
where from Eq. (|70|) , 
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^QCD 



with 

A^CD = A^exp(-(47r)V^o/(A)). (76) 

Solving the renormahzation group equations and substituting for g yields the following 
expressions for the mass and Coulomb renormahzation constants 

Zk{A) = Zk{Ai) I , "T" I =Zk{Ai)4^, (77) 



QCD ■ 



and 



47r r J A? 




(78) 



Lastly, we examine the effective renormalized potential between static color sources. This 
may be defined via Eqs. ([10|) and (^TD as 

V{k) = ^(^^^^(^^ = 4:TTa^s{k) .^g. 

rh K 

It is clear that it is the combination Zx(A)(7^(A) which is responsible for making a^s A- 
independent. For large A and Ai one obtains 

(.. (. , ZK{A)g\A) «cff(Ai) .sf.^ 

acff(A) = ttcflf(Ai)i^vTx^7T-^ = r- (80) 

Notice the power in the denominator which is present due to the rainbow-ladder nonper- 
turbative structure of the VEV of the Coulomb operator. Expanding Eq. (^Dj) permits a 
comparison to perturbation theory: 

aefr(A) = aeff (AO [l + ^/9o^^ log ^ + 0{g')\ . (81) 

In perturbation theory (with no light quarks) the coefficients in front of g^{Ai) should 
be equal to ANc — Nc/3 = 11 rather than 3/2/5o = 4iVc. The difference comes from 
the perturbative contribution due to emission and absorption of a transverse gluon, which 
involves iterating the Coulomb-transverse gluon vertex from He twice. This contribution is 
not present when one takes the expectation value of the Hamiltonian as done here. However, 
as stated earlier, perturbative contributions from propagating transverse gluons may be 
included, for example, in bound state perturbation theory and can be systematically included 
in our approach when the Hamiltonian is diagonalized in the quasiparticle basis. It should 
also be noted that such differences are of a screening nature, and thus are not expected 
to spoil the confinement mechanisms coming from summing the Coulomb-transverse gluon 
interactions. 
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D. Approximate Analytical Solution 

In this subsection we present an approximate analytical solution to the truncated renor- 
malized coupled gap equations for d{k), f{k) and uj{k)] Eqs. (|53[), (p^D, (p7|), respectively. 
The approximate solution is obtained by simplifying the angular part of the integrals over 
3-momenta. In each case the angular dependence is approximated by 

(k - q)2 ^ e{e - g2)fc2 + ^(^2 _ ^2)^2_ (g2) 

Next we assume that the renormalized solution of the gap equation can be written in the 
form 

uj{k) = 6{k — mg)k + 6{mg — k)mg. (83) 

Thus we assume that the gap function saturates to a nonzero value at low momentum. Once 
the FP operator d{k) and the Coulomb kernel f{k) have been obtained, the gap equation 
may be solved for uj{k) and the consistency of the ansatz for u may be checked. 

With the aid of these approximations the equation for the running coupling can be 
converted into differential form 



d'jk) [ (^ (ilf - l^d'ik)) for k < rUg 



(84) 



For k < irig the solution is given by 

d{kf d{fif 



which is well approximated by 



4 5 



(85) 



d{k) = ^ r- (86) 

This equation is trivially /i- independent. For large momenta, k > mg we approximate 
Eq. ( p^ by neglecting the terms of 0{m? Jk'^). In this case the solution is given by 

d{k) = ^^^ -r, (87) 

which also is /x- independent. Even though this solution is valid for k » irtg it may be 

matched continuously with the solution for kg < irtg if one chooses /i = rUg. The freedom in 

the renormalization of d{k) is now related to the choice of the value of d{k) at k = fi = rrig. 

It follows from Eqs. ( ^Sf ) and (|86D that there is a critical value of d{mg) = dc = 



47r-i/3/5/3o ~ 3.4414 for iVc = 3 which leads to d(k) oc l/yk for small k. Furthermore, 
this is the strongest possible IR enhancement admitted by the approximate solution. The 
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solution for d approaches a finite value for all other values of d{mg) less than dc- We shall 
see that this general behavior remains true for the full numerical solution as well. 

The corresponding solution for the function f{k) follows from Eq. (p3|). For k < rUg 

(with /i = rrig) one gets 



fik) 



fi^g 



1 + 



5/3o 
3{47r): 



and ioT k > m 



'91 



f{k) 



^d\mg) (^) 



fi^g) 



f{mg 



l+(l^rf^K)M$ 



f{mg 



, d{k) 
d{mg] 



, d{k) 
d{mg] 



(88) 



(89) 



The freedom in choosing the normalization for Zk{A.) is now reflected in the unspecified 
normalization constant f{mg). The maximal infrared enhancement of the Coulomb kernel is 
given by A;"'^/^ (/(^) oc l/A;^/^) if the approximate solution of Eq. ( P^D is used, or is given by 
^-15/4 (^f(^j^^ oc 1/k^^^) if the full solution for d(k), Eq. (^) is used. We note that a linearly 
rising Coulomb potential requires f{k) oc 1/k for small k. The exact numerical behavior of 
/ will be discussed in the next subsection. Lastly, if one substitutes the ansatz solution for 
the gap function Eq. (|83D into the gap equation (|67|) , one finds that it is indeed a solution 
up to terms of order 0{k/mg) for k < nig or 0{mg/k) for k > nig. 

To summarize, the approximate analytical solution leads to a running coupling (FP 
operator), d{k) which falls off logarithmically at large momenta and is enhanced at small 
momenta. The approximate solution indicates that there is only one critical value of the 
coupling for which the enhancement is maximal and given by d{k) oc 1/yk. This may be 
an artifact of the truncation of the series of coupled self-consistent equations. One expects; 
however, that the critical behavior is universal, i.e. near the critical coupling higher order 
corrections to the vertices in the Coulomb operator become irrelevant. 

The full Coulomb kernel becomes logarithmically suppressed at large momenta as ex- 
pected from an all-order resummation of leading logs. At the critical point and for low 
momenta it becomes enhanced over the perturbative 1/fc^ behavior and scales as k~^^/'^. We 
have thus obtained a tantalizing glimpse of the possibility of constructing a phenomenolog- 
ically viable truncation of QCD. 



E. Numerical Solution 

Encouraged by the near-appearance of linear confinement in the approximate analytical 
solution we proceed to a full numerical solution to the truncated renormalized coupled gap 
equations. The solution is obtained by mapping the gap equations onto a set of discrete 
nonlinear equations by placing all functions on a momentum space grid. We have found that 
numerical stability is enhanced if the grid is chosen carefully, in particular by preferentially 
populating the low and high momenta regions. The discrete gap equations were then solved 
with two independent solution algorithms. Both methods used an iterative procedure to 
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cycle through the three equations. Convergence was typically achieved in only a few passes 
since the analytical starting point of the last section is quite accurate. 

The numerical and approximate analytical solutions for the FP operator are shown in 
Fig. 4 for three separate values of dinig). This and subsequent figures are plotted in units 
of TUg which after renormalization is the only dimensionful parameter. Its value can only 
be determined upon comparison to a physical observable. It is clear that the analytical 
solutions are very accurate. Furthermore, the existence of a critical coupling appears to be 
numerically confirmed, with a value very near dc = 3.5. 



^10 




FIG. 4. Solution for the expectation value of the FP operator d{k). The two lower dashed lines 
correspond to analytical, approximate solution with d{mg) = 2.5 (lower) and d{mg) = 3 (higher). 
Boxes correspond to a full numerical solution for the corresponding values of d{mg). The numerical 
solution close to the critical point is shown by the open circles. The solid line corresponds to a fit 
to this numerical solution using the formula in Eq. (|^ 



The numerical solution near the critical point has been fit to the formula 

dc[^) ioT k < irin. 



d{k) 



dr. 



log(l+bd) 



for k > rrin 



(90) 



The fit yields dc = 3.5, a^ = 0.48, bd = 1.41 and q = 0.4 verifying the accuracy of the 
approximate analytical solution. Fig. 5 shows the Coulomb kernel function f{k)/ f{mg) for 
d{mg) = 2.5 and 3.0. Again, for d{mg) < dc the solution saturates at low momentum and the 
analytical approximation is quite accurate. The solution at the critical point is compared 

with 

' [jf-^ ^ for /c < TUg, 

^2iil±h^Y' ioTk>mg. 



f{k)/f{mg 



(91) 



log(fc2/m2+6j) 

The fit yields, a/ = 0.97, hf = 0.82 and Cf = 0.62. The low momentum behavior is found 
to be more enhanced than in the approximate analytical solution. The two fits to the 
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numerical solutions for d and / result in the following expression for the Coulomb kernel 

V{k) = f{k)d\k)/e: 



k^V{k)/f{mg 



(3.50)2 (!^) 



1.93 



for k < m 



3' 



8.07 log""-''"(A;7m2 + 1.41) log-"-^"(A;7m2 + 0.82) for k > m 



(92) 



10' 



10^ 



10 



10 



10 




FIG. 5. f{k)/f{mg). Curves as in Fig. 4. The numerical solution near the critical point (open 
circles) is fit to the formula given by Eq. (pT|) (solid line). 



At low momenta the effective coupling acs{k) (defined through Eq. (|75D) behaves very 
nearly as l/k"^. The fact that the power is not exactly —2 may be due to discretization error 
(a finer momentum grid does indeed bring the coefficient closer to —2) or the truncations 
employed in deriving the gap equations. In any event, as will be shown later the difference 
(roughly 3.5%) is completely negligible with regards to phenomenology. 

Assuming linear confinement {2ad + a/ = 2) gives [0] 



(93) 

530 MeV. 

630 MeV. 

These estimates of the scale are in accord with lattice computations of the adiabatic hybrid 
surfaces (thus is discussed further in Sec. IV.A) and with old glueball phenomenology |2^ 



6nb = {3.5f ml. 

Inserting the quark model value for the string tension, b = 0.18 GeV^ yields rrig 
Alternatively, lattice string tensions are typically 0.26 GeV^ [^, giving rrig ? 



The numerical and ansatz solutions for the gap function are shown in Fig. 6. We note 
the remarkable accuracy of the simple ansatz for u, the main difference being the smooth 
transition through the intermediate momentum region. Notice also that u approaches k 
very rapidly for large momentum. 
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FIG. 6. Comparison of the analytical approximation to uj{k) (dashed line) and the full nu- 
merical solution (solid line). 

Finally, the numerical stability of the solutions have been tested by varying the number 
of grid points. Of course this also tests the de facto numerical cutoff dependence of the 
results. The results are shown in Fig. 7. We find that the numerical results are stable to 
within a percent. Notice that this also confirms that all UV divergences have been properly 
subtracted. 
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FIG. 7. Numerical solution for d{k) near the critical point for 192 (circles), 288 (triangles) and 
384 (diamonds) grid points. 

F. Higher Order Terms 

We now address the issue of the neglected terms in the coupled gap equations. These 
arise, for example, from truncation of the rainbow-ladder sums, higher order corrections to 
the Coulomb vacuum energy, and from the terms generated by the Faddeev-Popov determi- 
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nant J . 



1. Vertex Corrections 



The truncation to the rainbow-ladder resummation for the Faddeev-Popov operator and 
the Coulomb kernel ignores higher order, 0[d^\ n > 2 corrections to the triple Coulomb- 
transverse-gluon vertex. Using the approximate analytical solutions for d{k) and uj{k) we 
estimate these contributions by evaluating the 0{(P) correction. From Eq. (^) it follows 
that the lowest order correction to the bare vertex is given by 

AF'^n. . V .^ - ^^ / ^^ [(k + l)-^T(q)r[k-^T(l)(k + l-q)]c^(k + l)d(k+l-q) ^, 
,)i,(k,q,k qj__y^_^ —^^ ________j . 

(94) 

We have evaluated this integral numerically and found that for all values of the external 
momenta the correction does not exceed a few percent. 

2. Second and Fourth Order Corrections to the Coulomb Kernel 

Recall that an operator product expansion for the Coulomb kernel has been defined in 
Sec. II. E, Eq. (0). We now employ the Swift equation (|33|) and the operator expansion of 
the Faddeev-Popov operator Eq. ( P3[ ) to derive an explicit expression for the terms in that 
expansion: 



K^^k, p; A) = 6{k + p)5"''ir(°) (k) + tgr^'' [k^^^^ (p)D(°) (k) + D(°) (p)K(°) (k)] : A'^(p + k) ■ p : 

n 



dg 



n n— 1 

5(p + k - 5] s,) : A^i (si) ■ p . . . A^"(s„) . (p - ^ s,) : (95) 

The term in the expansion of K which contains n gluons is weighted by a product of n — 1 
factors of D^^' and a single factor of K^^' . The additional contributions to the VEV of the 
Hamiltonian discussed in Sec. II. E, Sq and Sq\ come from terms with a product of n = 2 
and n = A normal ordered gluon fields : A"' : respectively. These are the only contributions 
which have a nonzero VEV after combining with the charge densities. The contribution to 
the gap equation is then obtained by taking the derivative of the VEV with respect to uj. 
As was discussed earher, an alternative method to derive the gap equation is to require that 
the off-diagonal (proportional to a^a^ or aa) portions of the one-body operators vanish. 
The second method would indicate that terms with n = 6, p : A^ : p, contribute to the 
gap equation as well since the four gluon fields from the two charge densities can contract 
with the fields from the kernel leading to an operator proportional to : A^ :. As discussed 
in Sec. II. E, the apparent difference in these two procedures is resolved if one notices that 
there are contributions to the gap equation which arise from the implicit dependence of the 

26 



Coulomb kernel on uo. In the second method, the contribution which would be associated 
with the n = 6 term in the operator product expansion of K is identical to the one from the 
derivative of the kernel in the n = A term contribution to the VEV. This was denoted 8(j 
in Sec. II. E. Similarly the term referred to as S\j in the discussion preceding Eq. (|54|) is 
identical to the contribution from the n = 2 term when the fields from the charge densities 
are contracted with each other. 

Adding all these pieces together yields, 



u'^{q) -ivifif 



d' + [4°^"(?)] + [4°^'^(^) + 4^"(^)] + [4^^^(^) + 4^"(^)] 



+[4'^'''(?)]-(? 



/.). 



(96) 



The four terms in the brackets are 0((i^), 0(<i^), 0{d^) and 0{d^) respectively; no other 
corrections exist. We test the importance of the higher order terms by computing the 0{d^) 
correction to the truncated gap equation. 




FIG. 8. 0{d'^) contributions to the gap function from £(^ 
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An example of a diagram contributing to £q is given in Fig. §. The explicit expressions 



m,K 



,(2),a; 



for £^^ and £^(^ are given below. 



S. 



{0),K 

c 



NJ 



iQ) = ^ 



dh dp 



8 J (27r)3(27r) 



uj{p) 



Tr[<5r(p)<5T(k)][(k + p)(5T(q)(k + p) 



/(k + p)d'^(k + p)d{k + p)d{k + p + q) + perm. 
(k + p)2(k + p)2(k + p + q)2 



(97) 



and 



'(2),^ 



S'c'^'il) 



dh dp 



Nl 

16 J (27r)3 (27r)3 



1- 



uj\q) 



[(q + k)(5T(p)5T(q)<5T(k)(q + p)] 



X 



uj{k)uj{p) 
/(q + k + p)d'^{q + k + p)(i(q + k)d{q + p) + perm. 



(q + k + p)2(q + k)2(q + p)2 



+ - 



N? 



dk dp 

(27r)3 (27r)3 



u{k) 



(k-q),5T(p)5T(k)5T(q)(p-k) 



27 



/(p — k)(i^(p — k)(i(p — k + q)(i(k — q) + perm. 
(p-k)2(p-k+q)2(k-q)2 



(98) 



Here the permutations refer to the other two ways of arranging the argument of f(f in 

Including these terms in the gap equations modifies the resuhs for d and / by strength- 
ening the IR enhancement somewhat. The result for the gap function is shown in Fig. |^. 
As expected, the change at higher momenta is minimal. However, we see that the (9((i^) 
terms do not modify u at low momentum either. This is because S'^''^ and S^'^ depend on 
the combination u}(j)i)/ui(j)2) — 1 which suppresses them in the IR limit (pi = P2). Our 
results are compared to lattice computations in Sec. IV. A. We stress that the result of Fig. 9 
should be considered preliminary because there are other 0{d'^) corrections (see the next 
subsection) that have not yet been included. 
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FIG. 9. Normalized instantaneous transverse gluon propagator, l/2w(A;). The dashed line is 
the solution to the leading order gap equation of Eq. (|67|); the solid line includes 0{d'^) corrections. 

The computation of the 0{d^) and 0{d^) corrections is progressively more difficult and 
is currently under investigation. These require the numerical solution of a self-consistent 
equation involving at least 8-dimensional integrals. However, since the 0{d'^) corrections 
are small we expect these higher order terms not to change the results significantly. 



3. Faddeev- Popov Contributions 

We now discuss the corrections due to the Faddeev-Popov determinants J'. We calculate 
the contribution to the gap equation for the determinant present in the kinetic part of the 
Hamiltonian, through the jr~^/^Hjr^/^ operators. This is given by 
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]- I di^J-^'^W{^)JW{^)J^'^ = ]-[ t/xn'^(x)n"(x) + Va. (99) 

Similarly, Vb is defined via the relation (see Eq. ^T]) 

Hc = \j d^dyp\x)Kab{^, y; A)p\y) + Vb. (100) 

A direct computation yields 

2 

Va = ^- j di^r'Tf5Tiy^),,{-xb\{y ■ D)-^|x,c) Vx, <5r(Vx).fc(xe|(V ■ /^)-'|x,/) V^^ 

2 

_ |_|rfxr^7-/5^(V,),,(x6|(V ■ /^)-^|x,/) v., 5r(V.),,(xe|(V ■ D)-i|x,c) V., . 

(101) 
or in momentum space 

+i / (0(0^(0^^"^^^' ^ + P)^"(P' ^ - k)/"'^/""MT(k)p] (102) 

(103) 

where D°'^{p, k) = (ap|y^|ck). We note that Va is similar to Christ and Lee's Vi; however, 
it is not identical because we have not Weyl ordered the Hamiltonian. 

Using the operator product expansion for the FP operator, these lead to terms propor- 
tional to : A^ : which add to the gap equation the following contribution 

2.^. N^ f dk dp rf(k + q)rf(k) ,2d{p + q)d{p) 



16 7 (27r)3(27r)3(k + q)2 k2 ' '^^'^' (p + q)^ p^ 
Nl f dk dp d{k + q)d{k) d{p + q) d{p) 



-[k6T{q)p][{k + q)6T{p + q + k)(p + q)]- 



16 J (27r)3(27r)3(k + q)2 k^ ' ^'^^^^'^^ ^^"^ ^^ ' ^ ' ^^^ ' ^^^ (p + q)2 p2 

The contribution of Va to the gap equation is IR-finite but UV-divergent thus will modify 
the gluon mass counterterm. A detailed numerical study of the full 0{d^) corrections to the 
gap equation will be presented elsewhere. 

IV. DISCUSSION 

As demonstrated in the previous section, the asymptotic behavior of the numerical solu- 
tion to the gap equations is V{k) ~ l/k'^, it thus appears that the methodology advocated 
in this paper is capable of describing quark confinement. The appearance of the confine- 
ment phenomenon hinges crucially on the choice of the variational vacuum which we use to 
construct the quasiparticle basis and on realizing that this choice also affects the interac- 
tion between these quasiparticles via the summed expression for the instantaneous Coulomb 
kernel. We now examine the implications of this success on confinement and the Gribov 
ambiguity. 
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A. The Confinement Potential 

The requirement that the gluon mass gap function uj{k) be cutoff independent gives rise 
to a mass scale which we call the gluon mass, rrig. The value of ca; at a particular momentum 
scale, say k = serves then as the underlying mass parameter of the theory. At the critical 
coupling dc the only free parameters in the gluon sector are Z^i^j) and the momentum scale 
itself, TUg = uj{0). Nonperturbative renormalization may be carried out by requiring that the 
Coulomb kernel reproduce the static QQ heavy quark potential as seen on the lattice (recall 
that Vqq is a renormalization group invariant quantity). In our approach this potential is 
given by 

HQCD\Qir/2),Qi-r/2)) = yQQ(r)|Q(r/2), Q(-r/2)), (105) 

In pure QCD, i.e. ignoring light flavors, the above eigenstate can be expanded in terms of 
multigluon states constructed from the quasiparticle operators acting on the \uj) vacuum. 
Schematically, 

|Q(r/2),Q(-r/2)) = Y.Ma^rb\r/2)S{-r/2)\u). (106) 

n=0 

where the quark creation operators refer to static sources. The Hamiltonian mixes states 
differing by gluon number; however, one expects that the mixing between such states to be 
suppressed by energy denominators due to the gluon mass gap, a;(0) 7^ (this is discussed in 
much more detail in Sec. V.A). This mass gap can be estimated from the difference between 
the lowest and excited adiabatic QQ potentials which have been calculated on the lattice 
T2\. One finds that this difference is AV{r ~ 1 fm) ~ 600 — 800 MeV. This is a natural 



estimate for uj{k) at low momenta. The implication is that the static ground state heavy 
quark potential may be accurately computed by ignoring extra gluonic excitations in the 
heavy quark system. (A calculation of the excited adiabatic potentials will be presented 
elsewhere.) Thus, to good accuracy, the static heavy quark potential V{r) is given by 

VQQ{r)--C^V{r)--C^J j^^e 1,2 = '^f J j^,e —^^. (107) 

It is useful to return to the approximate analytical solutions to the truncated renor- 
malized gap equations of Sec. III.D to illustrate how the different parameters enter. We 
have seen that at the critical point dc the solutions for the Faddeev-Popov operator and the 
Coulomb kernel are 

d{mg) (^j^ ioT k < TUg 
d{k) = { ,i^\) ior k>mg (108) 



1+3 ijd. 



and 



f(k)^{ /c'-.) fo,,>,„, (109) 



3,„„ fc2 y 



flog- 

■nig- 
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Since d{mg) = dc ^ 3.4414 is fixed, the potential has only two free parameters, the overall 
strength determined by firrig) and the mass scale set by rUg = u{0). These may be deter- 
mined by comparing with lattice computations of the Wilson loop. One finds firrig) ~ 1.0 
and nLg k, LS/rg. Here tq is the Sommer parameter of lattice gauge theory which is deter- 
mined to be roughly 1/430 MeV"^ Thus rUg k. 720 MeV. 

The same procedure may be followed for the numerical solution to the gap equation. 
Good agreement with the lattice static potential is obtained by choosing f{mg) = 1.41 and 
rrig = 1.4/ro = 600 MeV. The minimum in parameter space is fairly broad, for example 
firrig) = 1.09 and rrig = 1.6/ro = 690 MeV provides nearly as good a description of Vqq. 
The resulting potential (after numerically Fourier transforming to configuration space) is 
presented in Fig. ^. One sees that the numerically obtained static quark potential provides 
a reasonable facsimile of the lattice potential. This somewhat surprising result provides a 
posteriori support for the methodology advocated here. 
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FIG. 10. Static QQ ground state potential. The solid line is the full numerical solution for 
f{mg) = 1.41 and nig = 1.4/ro = 600 MeV. Data are taken from Ref. [22]. 

B. The Gribov-Zwanziger Horizon and the Gap Function 

As mentioned in Sec. II. B, the Gribov problem may be resolved by selecting a single 
gauge copy from the ensemble of Gribov copies by imposing the horizon condition of Eq. 
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(P3|). Furthermore, Zwanziger [0] has shown that the restriction to the fundamental modular 
region imposed by the horizon term is equivalent to having a low-momentum enhancement in 
the Faddeev-Popov operator over the perturbative l/fc^ behavior. This enhancement takes 
the form 

A;,A;,(E(0)i, - S(A;)i,), (110) 



where S is a regular function at the origin. Such a behavior is clearly an indication of 
confinement ^ , since the FP operator determines the static potential between color sources 
(cf. Eq. (0) or Eq. (^)). Comparing Eq. (|110|) with Eq. (|38|) shows that this is equivalent 
to the statement that d{k) is singular at the origin, 

lim l/d(fc) -> 0. (Ill) 

fc— »o 

The behavior of d at small momenta depends on one integration constant, d{fi). As we have 
shown earlier, d{k) approaches a finite value as A; ^ except when d{fi) = d^ where the 
Gribov-Zwanziger singularity develops. 

It is possible that the saturation of d{k) to finite values when d{^) < dc is an artefact of 
the rainbow-ladder approximation - we leave this as a matter of future investigation. For 
the present we simply require that the theory give rise to an enhancement of the FP operator 
at small momentum, this boundary condition then selects the coupling d{fi) = dc- 

In Ref. [0 the enhancement in the FP operator was obtained by adding the horizon term 
to the Hamiltonian via a Lagrange multiplier. The VEV of the new Hamiltonian was then 
computed in the bare vacuum, i.e with uj{k) = k; however, the horizon terms adds a mass 
term through an effective A^ operator whose strength is determined by the expectation value 
of the FP operator. This term has the effect of enhancing d{k) for small k. The equivalent 
of Eq. (|^) was then solved and the Coulomb operator was approximated by the square of 
the FP operator. 

In our approach the horizon condition was used to justify the expansion of the FP 
operator in a power series in gA. The resulting expressions for the FP operator were summed 
in the presence of a nontrivial mean field background (the variational vacuum) with the aid 
of the rainbow approximation producing IR enhanced FP and Coulomb VEVs. The success 
of this procedure demonstrates that the explicit effects of the horizon term may be ignored 
if one is willing to develop the quasiparticle spectrum and interaction self consistently. 

To further test this mechanism for realizing the Gribov-Zwanziger confinement scenario 
we compare our result for the gap function to that computed by Cucchieri and Zwanziger 
in SU(2) lattice gauge theory p^. In that paper, the authors measure the transverse and 



instantaneous gluon propagators in 'minimal Coulomb gauge'. They compare the numerical 
results to a functional form proposed by Gribov [|]: 



We call the scale appearing in this relationship, the Gribov mass, Mq- Cucchieri and 
Zwanziger found the the computed instantaneous transverse propagator agreed very well 
with this functional form but does not reproduce the normalization. 
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As discussed earlier, in our approach the transverse gluon propagator is suppressed due 
to an infrared singularity in the one-body gluon operator. Explicitly, the one body operator 
in the quasiparticle basis is given by 

r dk 
Hone-body = J2 J 7^-T^^(^)"^(k, A, c)a(k. A, c) (113) 



A,c 



with 

E{k) = uj{k) 



\^N^f dq /(k-q)rf2(k-q)l + (k.q)2 



4 J (27r)3 (k-q)2 uj{q) 



(114) 



The low-momentum enhancement of the kernel makes the integral infrared singular. Thus, 
as expected, gluons do not propagate. 

We note that the equal time transverse gluon propagator is not determined by E{k) but 
by uj{k) 

D'^ix) = lim(c.|T[A'^(x,t)A^(O,0)]|c.) = ^"^ (2^)3 t(fc) '^''^- ^^^^^ 

We have seen that the gap function obtained in Sec. Ill is rather flat at small momenta, 
even when some d^ corrections are incorporated into the gap equation. This is inconsistent 
with the lattice calculation for the same quantity, but as shown above, not inconsistent with 
the Gribov confinement scenario. The disagreement with lattice may be due to the use of 
the rainbow-lattice approximation and is being investigated. It is worth noting however that 
the effective gluon mass found by comparison to the QQ potential (or alternatively the D^'^ 
gluon propagator) is consistent with that found in the Coulomb gauge lattice calculations. 



V. IMPLICATIONS FOR THE CONSTITUENT QUARK MODEL AND 

PHENOMENOLOGY 

We now turn to an examination of the implications of the results presented here on the 
phenomenology of hadrons. Since this depends crucially on the explicit definition of hadronic 
states, we begin by searching for an efficient way to construct hadrons by specifying a new 
constituent quark model of QCD. The phenomenology of confinement is then analysed in 
light of the results of the last two sections. We conclude with a clarification of several 
open issues in the old constituent quark model and present a justification for the surprising 
efficacy of the quark model for light hadrons. 

A. Constructing Hadrons 

It is clear that constructing hadrons from the basis of free quasiparticles is futile if it is 
done perturbatively. A simple and natural way to avoid this pitfall is to choose a convenient 
form of Ho and diagonalise it nonperturbatively to obtain a basis of color singlet bound 
states. Bound state perturbation theory may then be employed to systematically include 
the effects of Hint- In our case the natural assignment for these operators is 

33 



Ho = J^^ {-la ■ V + /5m)^ + iy^xn^ - ^Jd^AV^A + ^ J d^dy p''{^)K^'\^ - y)p"(y) 

(116) 

and 

Hint = i/rfx[B2 + AV'a] -gJij^a--Aij + VA + VB + 

+ i/dxrfyp'^(x) [is:"''(x-y;A)-(5'^^irW(x-y)]p^(y) (117) 

The general philosophy is clear - Hq generates hadronic bound states; Hint incorpo- 
rates the corrections to these states due to transverse gluon exchange, three and four gluon 
interactions, and higher order contributions from the FP determinant and instantaneous 
confinement potential. It is worth stressing that Hq is still a field theory and hence is con- 
siderably tougher to solve than old fashioned quantum mechanical quark models. But there 
are substantial advantages to adopting this approach. Foremost is that Hq + Hint is QCD. 
Furthermore, Hq is relativistic and incorporates gluonic degrees of freedom. Thus it is pos- 
sible to examine glueballs, hybrids, and other gluonic phenomena in a coherent fashion. The 
utility of the rearrangement made in Eq. ( [L17| ) lies in the use of the variational vacuum to 
construct a phenomenologically viable basis of quasiparticles. This has the direct effect of 
greatly improving the Fock space convergence of any observable. As we have seen, it also 
automatically generates the correct static potential upon which to construct hadrons. We 
have previously mentioned that Ho generates states which are infrared divergent if they are 
not color singlets (hence these are removed from the spectrum). Conversely, all color singlets 
are IR finite. Thus the basis generated by Hq contains no spurious color nonsinglet states 
which would have to be removed by laborious iteration of Hint and, in fact, is expected 
to provide a reasonably accurate starting point for hadronic spectrum computations. As a 
practical note, the physics of the variational vacuum may be accurately approximated by 
simply using dressed quarks and gluons when constructing hadrons. The constituent masses 
are roughly 200 MeV and 600-800 MeV respectively. Finally, the spectrum generated by Hq 
is spin averaged in the sense that it only incorporates spin effects from relativistic corrections 
to the Coulomb potential. Full spin splittings come from Hint- 

An important implication of this approach is that the rapid convergence of the con- 
stituent quark model Fock space expansion has a natural and simple explanation. All of the 
corrections induced by Hint (for nonexotic states) involve the transfer of a virtual transverse 
gluon. Since these are quasiparticles in the variational vacuum, the relevant perturbative 
diagrams are suppressed by the mass gap between the regular and hybrid states. This sim- 
ple feature of QCD in Coulomb gauge has important phenomenological consequences. For 
example, it implies that the Fock space expansion converges quickly because state mixing 
involves the creation of massive gluonic (or quark) quasiparticles. Recently, lattice data 
has appeared which confirms this picture. Duncan et al. [^ have constructed a simple 
relativistic quark model of B mesons by considering a light relativistic quark (with kinetic 
energy Vfc^ -|- m?) moving in the lattice T potential. Detailed comparison with lattice B 
data demonstrated the high accuracy of the model. The point which is relevant for our 
discussion is that the lattice T interaction (recall that this is equivalent in principle and 
in practice to K'^'^^) should receive corrections due to the light quark when applied to B 
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mesons; see Fig. ^. The fact that these corrections are not important demonstrates that 
they are suppressed, in agreement with the above arguments. 




FIG. 11. The leading light quark correction to the confinement potential 

B. Confinement in the Constituent Quark Model 

One of the benefits of Coulomb gauge is that it makes the source of confinement clear: in 
the heavy quark limit quarks and transverse gluons decouple and the quark-antiquark quark 
interaction arises solely from the instantaneous Coulomb operator. This rigorous result has 
several significant implications for hadronic phenomenology. 

First we simplify the situation by noting that higher order terms such as shown in 



Fig. |I2| are suppressed due to the arguments espoused in the previous subsection. Thus, 
the dominate interaction between static color sources is the leading kernel in the Coulomb 
interaction, K^^\ As we have seen, this kernel is essentially identical to the lattice Wilson 
loop result, so this conclusion is supported a posteriori. 
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FIG. 12. Examples of Higher Order Corrections to the Heavy Quark- Antiquark Interaction. 

This simple statement carries wide repercussions. For example, a longstanding corner- 
stone of quark model phenomenology is that confinement is 'scalar'. What this means is 
that the interaction between quarks is assumed to be 
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^y'V^V^(x)ir(x-y)^^(y). (118) 

This form (as opposed to 'vector' confinement ip'^ipKip'^il)) is supported by a comparison of the 
predicted spin sphttings in heavy quarkonia with data [^. However, the resuhs presented 
here make it clear that this conclusion is naive. The interactions between color sources is 
more complicated than the simple facsimile given in Eq. ( p.l(j| ). As we have seen, the leading 
interaction between quarks is given by i^*-'^-* - and this has the form of vector confinement. 
What is taken as evidence of the scalar nature of confinement is in fact quarkonium spin 
splittings which are generated by nonperturbative mixing with intermediate hybrid states 
via Hint- That this more complicated (and correct) picture may look 'scalar' has been shown 
in Ref. m 

Another simple conclusion of the picture being developed here is that the confinement 
potential between color sources scales as the quadratic Casimir. This follows from the 
observation that the dominant contribution to the confinement potential is given by the 
leading kernel and that the color structure of this kernel is i^^j, = ^afe-R'*-"^. The fact that 
Casimir scaling of the Wilson loop potential has been observed repeatedly |2^ may be 
taken as a successful prediction of our methodology or may be used as further proof that 
the diagrams of Fig. |T2| are suppressed with respect to K^^'^ . 

The methodology presented here allows for the resolution of several open, but often 
ignored, ambiguities in the constituent quark model. For example, it is often stated that 
the linear potential is built from the exchange of infinitely many gluons. One may then 
ask why the one gluon exchange potential is retained as an important part of quark model 
phenomenology. Indeed the split between one gluon exchange color Coulomb and hyperfine 
forces and the multiexchange linear force is necessarily ambiguous. The resolution to this 
issue is transparent in Coulomb gauge: 'one gluon exchange' is part of Hint and is due to 
noninstantaneous transverse gluon exchange. The instantaneous central portion of the quark 
model should consist of a linear term in addition to the running resummed 'Coulomb' term 
of Eq. (|9^) . No ambiguity exists because of the separation of instantaneous and transverse 
degrees of freedom inherent in Coulomb gauge. 

Another problem with the old-fashioned CQM has to do with the previously mentioned 
assumed scalar nature of confinement. Unfortunately, scalar confinement implies that if 
mesons are bound by a linear potential, baryons are antibound |2^. This is clearly an 



intolerable situation which is routinely ignored by CQM practitioners. As we have seen, 
the resolution is that confinement acts as the time component of a vector rather than as a 
scalar, and no inconsistency exists between mesons and baryons. 

C. Constituent Gluons and Strong Decays 

We illustrate the power of our approach by considering the vexing problem of strong 
decays in hadronic physics. 

The strong decays of hadrons has been, and remains, a mystery of soft QCD. The naive 
perturbative assumption that the decay proceeds via one gluon dissociation (Fig. |l3|b) is 
proven incorrect by direct comparison with experiment [^. The only reasonably successful 



phenomenology is provided by the '^Pq' model [0, where quark pairs are assumed to appear 
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with vacuum quantum numbers over all space. This is clearly an unacceptable situation, 
especially given the ubiquity of hadronic decays and the fact that they provide a window 
into the dynamics of glue at low energy. 

We now examine the predictions of the new quark model presented here. To lowest order 
in KqcD/fT^g and to all orders in the coupling, the only diagrams which contribute to meson 
decay (here all mesons are assumed to have Fock expansions which are dominated by the 



quasiquark-quasiantiquark components) are shown in Fig. |I3|. The left figure is contained 
within Hq and is therefore the leading diagram. The central and right figures contribute at 
0{KQCD/'^g) and are generated by Hint- 

Diagrams (a) and (b) with perturbative gluons or model potentials in the intermediate 
states have been previously examined as possible sources of hadronic decays in Ref . |^ . The 



authors noted that diagram (a) is strongly suppressed with respect to (b) due to momentum 
routing through the pair production vertex (this diagram is zero in the nonrelativistic limit 
when a delta function potential is in place. It is strongly suppressed with a 1/q^ potential). 
The other class of diagrams considered in Ref. |^ was that generated by a phenomenological 
scalar interaction given in terms of scalar confinement (cf. Eq. (|118|) ). This is, of course, an 
ad hoc microscopic realization of the ^Pq model. What was found was that this diagram (like 
diagram (a) but with scalar as opposed to a vector vertices) was much larger than diagram 
(b). 
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FIG. 13. Leading Order in 1/mg Meson Decay Diagrams 

These conclusions imply that the ^Pq model would emerge in a natural way from our 
methodology if diagram (c) produced light quark pairs with scalar quantum numbers. Di- 
agram (c) is generated by the product of K^^^ and Hgg terms in Hint (see Eqs. (pO|) and 



^)) and is roughly given by A ■ 'Vip'^a ■ Aip. Once the vector potentials are contracted (or 
better yet, the sum over intermediate hybrid states is made), the resulting operator is of the 
form ilj'^[cr6T'V]tp, very nearly equal to the long-assumed ^Pq vertex. Thus we have obtained 
a viable microscopic description of hadronic decays. The implications of these observations 
will be explored in a future publication. 

D. Light Quarks and the Constituent Quark Model 

The utility of the CQM for heavy quarkonium is not in doubt. However, its appar- 
ently successful extension to light quark states is unexpected and surprising. We seek to 
understand this observation in this subsection. 
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The major feature of light quark physics is spontaneous chiral symmetry breaking. One 
may regard this as occurring due to the appearance of a quark-antiquark vacuum condensate. 
The interactions which generate the condensate are typically associated with an effective 
instanton interaction [^ or the confinement potential 0-0. (In our approach the driving 
kernel would be K^^'). Regardless of the particular mechanism which causes attraction in 
the scalar channel, a massive constituent quark is the necessary outcome. Indeed, while bare 
quarks may become very light or massless, the relevant quasiparticles saturate at roughly 
200 MeV as the bare quark mass is reduced [0,|T^]. This, at least partly, explains the 
apparent success of the nonrelativistic portion of the CQM. The agreement is also enhanced 
by the empirical accident that the expectation value of y/p"^ + m^ is very close to p^/2mcQM 
in typical hadronic states. More important than this; however, is the nature of the central 
potential itself when the bare quarks are light. As discussed above, effects due to one gluon 
exchange are suppressed by powers oi Kqcol^ng. Thus the main effect due to light quarks is 
the presence of intermediate quark loops in the instantaneous interaction (Fig. |14[). These 
diagrams cause string breaking which is an important feature of QCD. However, as Isgur has 
argued [Q , the main consequence of this is simply to renormalize the string tension. Thus 
light quark loops have little effect on the phenomenology arising from Hq. The conclusion 
is that structure of the new CQM which we have laid out is essentially unchanged for 
light quarks. Furthermore, even the simple nonrelativistic approximation may retain some 
validity for massless bare quarks. 




FIG. 14. Light Quark Loop Correction to the Heavy Quark- Antiquark Interaction. 



An explicit demonstration of how the CQM may emerge was given in Ref. |T^. This 
paper assumed a simple contact interaction in place of the full Coulomb kernel. Stan- 
dard many-body techniques were used to obtain chiral symmetry breaking and constituent 
(quasiparticle) quarks. It was then demonstrated that the vector meson - pseudoscalar me- 
son mass splitting follows a form essentially identical to that of the CQM hyperfine splitting 
when considered as a function of the constituent mass. Nevertheless, the mass splitting was 
clearly driven by chiral symmetry breaking when considered as a function of the current 
quark mass, thereby demonstrating that the pion may be viewed as both a pseudoGoldstone 
boson and as a quark-antiquark bound state. The new quark model presented here provides 
an explicit microscopic realization of the contact model employed in Ref. ||T^ and it will be 
of interest to verify the findings of that work. 
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VI. CONCLUSIONS 

In the paper we propose a new way to organize QCD which is appropriate for low energy 
hadronic physics. The starting point is chosen to be the QCD Haniiltonian in Coulomb 
gauge because this gauge is most directly apphcable to bound state physics - the degrees 
of freedom are physical and an instantaneous potential existsQ. The instantaneous Coulomb 
potential may be incorporated into Hq (as is done in atomic physics) and a viable bound 
state perturbation theory may be constructed. This simple step already obviates one of 
the severe problems of perturbative QCD in describing hadronic properties, namely that of 
ill-defined asymptotic states. 

While the division of the QCD Hamiltonian is a simple task, it is essentially meaningless 
because the degrees of freedom represented in Hq are partonic. Thus building bound states 
would be a frustrating exploration of the depths of Fock space rather than the preliminary 
step for bound state perturbation theory we desire it to be. The experience provided by the 
constituent quark model points the way out of the impasse: appropriate (constituent) degrees 
of freedom must be employed. The problem in the past (cf. constituent quark models, bag 
models, flux tube models, etc) has been in finding a way to introduce effective degrees of 
freedom in such a way that the connection to QCD is not destroyed. Herein we present 
one way to do this which is based on experience gleaned from many-body physics often 
used in phenomenological models e.g. the Nambu-Jona-Lasinio model [Q. Specifically, a 
canonical transformation to a quasiparticle basis which is defined with respect to a nontrivial 
variational vacuum is made. The theory remains QCD but is given in terms of a useful and 
tractable basis. Although the vacuum state is necessarily an ansatz, this does not vitiate 
the construction - in principle any basis may be used, we merely seek an efficient one, and 
the vacuum itself may be systematically improved with standard techniques. 

One finds a welcome complication when these ideas are applied to nonabelian gauge 
theory: the interaction which is needed to define the vacuum ansatz and the quasiparticle 
spectrum (via the gap equation) itself depends on the vacuum. Thus the fundamental 
quasiparticle interaction and the quasiparticles themselves are inextricably interdependent. 
Solving the gap equations requires the evaluation of the specific functional dependence of 
the quasiparticle interaction on the vector potential. We have chosen to do this within the 
rainbow ladder approximation. There are several important points to make at this stage: 
(1) the rainbow ladder approximation may be improved at will, (2) the approximation is 
accurate in the large Nc limit, (3) the approximation is accurate in the infrared limit, (4) 
the approximation is justified a posteriori. Lastly, although the approximation cannot yield 
nonperturbative results, true nonperturbative physics may be generated when the resummed 
kernel is incorporated in the nonlinear coupled gap equations. Doing so reveals a pleasant 
surprise: the emergence of the confinement phenomenon. 

While it is gratifying that color confinement is produced by our approach, this result 
would be useless if it did not match phenomenology. The fact that the effective quasiparticle 



^We note that it is also useful for QCD at finite temperature because a special frame is automati- 
cally selected and because counting degrees of freedom is an important aspect of thermodynamics. 
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potential matches the lattice static quark potential very well points to the general utility of 
our method. Thus Eqs. ( |116|) and ( |117|) represent much more than a simple reordering of the 
QCD Hamiltonian. By building Hq as an effective Hamiltonian describing the interactions 
of quasiparticles on a nontrivial vacuum we are able to establish contact with the constituent 
quark model and derive confinement. That both of these emerge in our formalism bodes well 
for the future success of Hq as a robust starting point for detailed hadronic computations. 

An important test of any new method in QCD is its ability to provide insight into a variety 
of phenomena. We have tried to demonstrate the robustness of our method in this regard. A 
vital aspect of this robustness is the emergence oi Aqcd/^^q as an expansion parameter. This 
provides the justification for gluonic Fock space truncation, for the validity of the leading 
static Coulomb kernel K^^\ and for the applicability of the static kernel to light quarks. 
Indeed, the method strongly hints as to why the constituent quark model works for light 
quarks. To summarize, quarks never become truly light (but saturate at constituent masses), 
the static kernel is not strongly affected by the presence of light quarks, and parameter 
freedom in the definition of the quark model allows for an accurate reproduction of the 
relativistic quark kinetic energy and the chirally-driven meson hyperfine splitting. 

The ideas we have presented have had a long period of development starting with Gribov's 
speculation that confinement may arise naturally when resolving the gauge copy problem. 



In the early 1980's Finger and Mandula [^, Adler and Davis 0, and Le Yaouanc et al. |TD 
all considered the generation of constituent quark masses and spontaneous chiral symmetry 
breaking with simple (often of the form given in Eq. [I18| ) models of QCD. The issue of 



renormalization was taken up by these papers and in Refs. [|T^,0,^ . 

The work which is closest to ours is that of Zwanziger [^ and Swift [^ . As discussed in 
Sec. IV.B, Zwanziger has shown that the imposition of the horizon condition implies that the 
Faddeev- Popov propagator is enhanced in the infrared. As we have stressed, an enhancement 
of the FP propagator is sufficient to cause confinement. In Ref. Zwanziger has shown that 
adding the horizon term to the Hamiltonian produces an effective gluon mass which in turn 
induces the desired enhancement of the VEV of the Faddeev-Popov operator. Zwanziger 
then makes several simplifying assumptions to arrive at an estimate for the Coulomb kernel. 
Chief among these are an assumed form for the gluon dispersion relation, a simplified version 
of the Faddeev-Popov propagator integral equation, and the approximation K ~ (P{k)/k'^. 
The end results are similar to ours; Zwanziger obtains d ~ k~^^^ (we get k^^^"^) and V ~ r^/^. 
Our analytical approximation gives V ~ r^/^ while the numerical solution is very nearly 
linear. 



The work of Swift E^ is very similar to ours in philosophy. In fact our self-consistent 



equations for the leading rainbow-ladder gap equations, which were derived in the Hamil- 
tonian formalism, agree with those of Ref. ETI], which were derived in the Greens function 



formalism. However, a difference occurs in the renormalization of the mass gap equation: 
we find that only one subtraction is necessary to render the equation finite. Thus no coun- 
terterm proportional to AV^A is required. This is due to the logarithmic suppression of 
the potential at large momenta. 



The main difference between the current paper and Ref. [ 21 1 is in the analysis of the gap 
equations. We have obtained very good analytic and full numerical solutions to the coupled 
gap equations. This was not attempted in Ref. Ell]; however, the author did examine 



the small momentum behaviour of the Faddeev-Popov propagator and the Coulomb kernel 
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assuming a particular form for the gap function. His FP and gap functions agree with our 
analytical estimates, however his solution for the potential has an unexpected imaginary 
portion. We believe this is due to an approximation which generated a confinement potential 
which was more singular than l/k"^. 



A preliminary exploration of the work presented here was undertaken in Ref. []13|. This 
reference neglected the FP determinant and higher order contributions to the gap equation. 
Furthermore, the full Coulomb kernel was drastically simplified by taking it to be the sum- 
mation of the one loop expression for K. Despite these simplifications and assumptions the 
resulting potential was similar to that obtained here. This is perhaps an indication of the 
power of the coupled quasiparticle/vacuum approach. 

We regard the present work as a promising start to the construction of a new quark model 
of the strong interactions in particular with regard to the treatment of gluonic degrees of 
freedom. Future projects include the evaluation of all d'^ correction terms to the gap equation 
which are needed to test the Gribov-Zwanziger gluonic quasiparticle spectral function. We 
also intend to evaluate a broad swath of the meson, baryon, glueball, and hybrid spectra, 
and to compute the heavy hybrid adiabatic energy surfaces. The latter two test the utility 
of the gluonic quasiparticles as effective degrees of freedom and will probe the structure of 
K^^^ and K^-^\ General considerations (and explicit lattice evidence) lead one to expect 
that glue behaves as a collective stringlike degree of freedom at large distances. We expect 
the gluonic quasiparticles to provide a useful description of glue at small (less than 1 fermi) 
and intermediate (1-2 fermi) scales. It will be interesting to see if the formalism presented 
here allows for effective stringlike behavior at large distances. Finally, rj — rj' mixing is 
a longstanding issue in soft QCD related to the anomaly, vacuum structure, the gluonic 
content of mesons, and instantons. Examining this problem should prove informative for 
the further development of the quark sector of our theory. 
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